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SUMMARY 


A family of FORTRAN IV subroutines was developed to calculate the 
thermodynamic and transport properties of eight fluids at pressures 
from 0.1 to 200 atmospheres (to 100 atm for helium) and from the triple 
point to 300 K for neon; to 500 K for carbon monoxide and oxygen; to 
600 K for methane and nitrogen; to 1000 K for argon and carbon dioxide; 
and from 3 to 500 K for helium. This paper presents some of the histor- 
ical basis for the calculations and the equation of state selected. An 
operations sheet with a detailed description of pertinent input param- 
eters and a discussion of the results is included in the text. A brief 
description of each subroutine is presented in appendix A. The master 
subroutine GASP calculates the remaining state variable when given any 
two of pressure* density, and temperature as input. In addition, any or 
all of the thermodynamic and transport properties - enthalpy, entropy., 
h specific heats (Cp, Cv) , sonic velocity, viscosity, thermal conductivity, 

surface tension - can be obtained. A call to GASP with pressure and 
,5 enthalpy or entropy as input will also generate the other properties. A 

special technique is provided to estimate the thermal conductivity near 
the thermodynamic critical point. 

The GASP package was developed to be used with heat transfer and 
fluid flow calculations. Because of this intent, the package was 
written as a family of independent subroutines with a separate subroutine 
for each specific property or set of related properties. Thus a user 
with limited needs can work with only those subroutines required for his 
program. 

Cryogenic fluids are used for cooling equipment, preservation, and 
destruction of biological specimens, heating or cooling other fluids, as 
modeling fluids, and in many cases as the primary test fluids in heat 
transfer and fluid dynamics research. Carbon dioxide is used in many 
experiments because of its "room temperature" critical temperature. 

Cryogenic fluids are used in biology and medicine for the preser- 
vation and destruction of tissue. Skin cancers and growths have been 
successfully treated by spraying liquid nitrogen on the defective 
tissue (1) , Thousands have been treated since 1962 for involuntary move- 
ment disorders (parkinsonism, dystonia musculorum deformaus, intention 
tumor, and tortioeollls) , using stereotoxic cryo thalamectomy (2). Gyne- 
cological disorders have been treated using cryogenic cannula techniques 
with a variety of freezing surface configurations (3) . 

A variety of techniques using cryogens are available for the pres- 
ervation of whole blood and its constituents for periods of years (4) . 
Cryogens have been used to preserve foodstuffs, such as fish, poultry, 
red meats, fruits and vegetables (mushrooms, tomato slices, peas, beans, 
berries, sliced peaches, etc.) (5 and 6). 



2 


Currently much effort is being devoted to whole organ preservation. 
The recent work of Lehr (7) indicates that the kidney may be able to sur- 
vive freezing when microwave energy is used to revitalize the organ. In 
cryoimmunology , destructive freezing has been shown to stimulate the pro- 
duction of antibodies (8) . 

Cryogenic fluids are used in the freeze drying of foods, micro- 
organisms, pharmaceuticals, and biologicals (9). Hydrogen and helium 
bubble chambers are used in nuclear physics (ref. 10). Ultrahigh vacuums 
can be achieved through cryopumping (ref. 11). Liquid oxygen (LOX) is 
used at hospitals; in increasing the capacity of sewage plants, and in 
1968, 70% of all LOX produced was consumed by the steel industry 
(ref. 12). High speed trains are proposed which require cryogenically 
cooled superconducting magnets (ref. 13). Superconducting power trans- 
mission equipment will require cryogens, helium (perhaps hydrogen), as 
well as cryogenic conditioning equipment. Ejectors operating in the 
cryogenic regime can substantially augment the capacity of refrigeration 
systems (ref. 14). Cryogens are used to produce a low background noise 
environment for infrared detection equipment (ref. 15). Cryoelectronics, 
superconducting motor-generators , and thermonuclear power are all depend- 
ent on cryogens . 

The problems associated with the liquefaction, storage, and gasifi- 
cation of LNG (liquid natural gas) and LPG (liquid petroleum gas) are 
currently being studied (ref. 16). 

Cryogens are used in metal forming cryoquenching , and for producing 
higher quality machine tools (refs. 17, 18). 

This list of applications is by no means complete and review of 
other cryogenic applications and several cited herein can be found in 
references 19, 20, and 21. The point is simply that each of these appli- 
cations requires not only a knowledge of the thermodynamic and transport 
properties of the fluid but an efficient readily accessible technique for 
evaluating these properties. 

GASP is a FORTRAN IV family of subroutines which was developed on a 
7094-7044 DCS machine (it is also available for the 360 system) to deter- 
mine the thermodynamic and transport properties of eight fluids (helium, 
methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide) 
and be used in conjunction with other research programs. The GASP prop- 
erty package is subdivided into subroutines and functions for computing 
the individual properties. The individual subroutines may be used inde- 
pendently of subroutine GASP if the user’s requirements are of a more 
restricted nature (i.e., memory or time limited). Appendix A is a de- 
scriptive breakdown of all the functions and subroutines used by the 
master subroutine GASP, which is listed in appendix B. 

The basis for the program is the work of Bender (ref. 22) and 
Strobridge (ref. 23). Strobridge curve fit the available nitrogen data 
with a modified Benedict-Webb-Rubin (BWR) (ref. 24) equation of state. 
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This work represented a major advance in establishing an equation of 
state for cryogens. The excellent hydrogen data of Roder and Goodwin 
(ref. 25) were fit using a form similar to that of reference 23. In the 
early days of work with cryogenic hydrogen such an equation represented 
a milestone to those working in research and development. 

Bender (ref. 22) added a new constraint to the problem of curve- 
fitting PVT data and required subcritical data to satisfy the Maxwell- 
Phase Rule. The use of this constraint, although not well defined in 
reference 22, can give saturation properties directly from the equation 
of state. Bender (ref. 22) fit selected PVT data, subject to the 
Maxwell-Phase Rule constraint, for five fluids: methane, nitrogen, 

oxygen, argon, and carbon dioxide from the triple point temperature at 
1/2 atm to approximately 600 K and 200 atm to within the accuracy of the 
experimental data. The form of Bender's equation (ref. 22) is based on 
that of Strobridge (ref. 23). Because ..the rules for obtaining saturation 
properties were not clear from referenq^ 23, the saturation properties of 
these five fluids were calculated herein from curve fits . 

fk 

Hust and Stewart (ref. 26) used the principle of corresponding 
states and the equation of Strobridge (ref. 23) to give thermodynamic 
property values of carbon monoxide. McCarty and Stewart (ref. 27) used 
a modified form of the Strobridge equation and corresponding states to 
fit the PVT data for neon. The neon data were subject to the usual 
critical constraints and one additional constraint; namely that the equa- 
tion must fit the data for the Joule-Thompson inversion locus. 

Mann (ref. 28) further modified Strobridge *s equation to fit the 
available helium data, much of which was compiled by Mann and Stewart 
(ref. 29). Since that time, R. D. McCarty (ref. 30) has obtained a 
more exacting fit to the data; furthermore McCarty's fit yields good 
values of S^p/ST^ for the fluid helium. This is not the case for the 
equations of Strobridge (ref. 23), Mann (ref. 28), Hust and Stewart 
(ref. 26), and McCarty and Stewart (ref. 27). Bender's equation 
(ref. 22) also gives good values of 3^P/3T^ for the five fluids fit 
therein. 

The equation of state used in this report is a modification of 
Bender's original equation (ref. 22) to include the modifications to the 
Strobridge equation (ref. 23) for the fluids, helium, neon, and carbon 
monoxide. The primary references used in GASP for thermodynamic and 
transport property calculations are presented in table I. The equation 
of state is presented in appendix A, under SUBROUTINE PRESS, and the 
coefficient matrix for all fluids is presented as table II. 

The Bender equation of state (ref. 22) has been compared to the 
P-p-T data of several investigators. Comparisons to data were made 
for methane, carbon-dioxide, oxygen, argon, and indirectly for nitrogen 
and helium. In figures 1 through 12, the relative errors in the calcu- 
lated P , p , or T are illustrated as functions of the reduced property 
for methane and oxygen. For example in figure 2, the relative errors in 
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pressure (P exp “ Pcalc)/Pexp are plotted as functions of the reduced 
pressure P exp /P crit . 

The essence of these figures and other comparisons can be partially 
summarized in table III. However, one should consult the figures in 
order to more fully understand the deviations between the state equation 
and the data. 

The results were generally in good agreement with the collection of 
data in Dinn (ref. 32). The fit of the helium data given by McCarty 
(ref. 30) indicates the helium results to be off as much as 4% in PVT and 
up to 20% in some derived properties in the near critical region, indi- 
cating the helium fit, used herein, will need further modifications. 

Even the neon densities are not in complete accord with reference 27, 
and the derived values deviate somewhat from those of reference 27, indi- 
cating that these data will also need further modifications. 

The GASP property package differs from the various techniques of 
references 22 through 30 in that it was developed to be used in fluid 
flow and heat transfer calculations. As such, there are independent 
calls for the three state variables pressure, density, or temperature 
(see table IV OPERATIONS SHEET). In addition, temperature and all the 
other properties can be obtained as a function of pressure and enthalpy, 
or pressure and entropy which are often used in forced convection 
studies. Selected examples of properties calculated using GASP are 
given as figures 14 and 15. 

While enthalpy and entropy are available in references 23 to 28 
the specific heats (Cp, Cv) , sonic velocity, viscosity, thermal conduc- 
tivity, and surface tension were not computed (good values of Cp and 
Cv are given for helium in ref. 30). One of the major reasons that 
the specific heats were not given is the significant deviations in 
3 2 P/9T 2 at near critical temperatures which generally lead to very 
large values of Cv. For example, using the Strobridge equation, the 
calculated values of Cv for nitrogen are illustrated in figure 13. 

After much effort in this area, the authors finally decided to use nu- 
merical techniques to determine AU/AT = Cv and AH/AT = Cp, for p/p c 
cutoff value that varies with each fluid. See figure 16 and SUBROUTINE 
CPPRL, appendix A. The viscosity and thermal conductivity computational 
methods were adapted from references 33 and 34 for the gaseous regime 
and from references 35, 36, 37, and 38 for the near-critical regime. 

The surface tension calculation was adapted from reference 39, and modi- 
fied to fit the date given in reference 40. A description of the equa- 
tions for these transport properties is included in appendix A. 

The symbols used herein are given in appendix C. 
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APPENDIX A 

FUNCTIONS AND SUBROUTINES CONTAINED IN SUBROUTINE GASP 

If the user wishes to disassemble GASP into its family of subrou- 
tines, the user should contact the authors at Lewis Research Center, 
Cleveland, Ohio 44135. 

In all cases , the units are specified by KU and the region by KR, 
see sketch, table IV. If KR is either specified as 1 or returned as 1 
the saturated liquid and vapor properties are computed. All saturated 
properties become available through the labeled common statement COMMON/ 
PROPTY/etc. Three of five choices of KU available to the user are 
given in table V. 

The property routine GASP will warn the user if the input or calcu- 
lated values are out of range, if an iteration has not converged, or a 
value has been extrapolated; however , calculation will continue. Calcu- 
lations will be terminated if the user fails to specify the name of the 
fluid, see SETUP (NAMGAS) . 

A flowchart of subroutine GASP and the routines called directly by 
GASP is included in appendix D. 


BLOCK DATA - A Data Subroutine 

This routine contains the message which lists the fluid used and the 
associated critical constants, e.g., 

THERMODYNAMIC AND TRANSPORT PROPERTIES FOR C02: 

PC - 72.869 ATM, TC = 304.21 K, ROC = .464 GM/CC 

Also included are the units conversion factors to convert the user’s 
units as specified by KU to those which are internally consistent. 

KU - 1 gives units which are the same as in the program. 

SUBROUTINE SETUP (NAMGAS) 

This routine uses NAMGAS, a 2 or 3 letter Hollerith code tabulated 
below, to select the coefficients for a particular fluid and overlay them 
into the common blocks needed by GASP. The coefficients include: 

(1) Equation of state coefficients 

(2) Saturated vapor law coefficients 
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(3) Specific heat at zero pressure coefficients 

(4) Transport coefficients 

(a) Viscosity y - y , y 

(b) Thermal conductivity X - X , X 

(c) Force constants 

(d) Surface tension a 

(5) Region delimiters on P, T, p, H 

(6) Region delimiters for switching to numerical Cv calculation 


Fluid 

Hollerith code 

Fluid 

Hollerith code 

Helium 

HE 

Nitrogen 

N2 

Methane 

CH4 

Oxygen 

02 

Neon 

NE 

Argon 

AR 

Carbon -monoxide 

CO 

Carbon-dioxide 

C02 


The coefficients for items 1 to 3 and 5 are listed in table II. The ref- 
erences for the other coefficients are listed in the fluid reference 
table I. 

SUBROUTINE CHECKS (A multiple entry subroutine) - 
ENTRY DCHECK(KU,D) 

At this entry point the density D is converted to 
grams/cnr and checked to see if D is out of range. 

ENTRY PCHECK(KU ,KR,P) 

This entry converts the pressure P to mega-Newtons/ 
(meter) ^ and checks for out of range. If KR=1, P is 
checked for out of saturation range. If P is out of 
range, the program writes an out of range note and 
continues . 

ENTRY TCHECK (KU , KR , T ) 

This entry converts the temperature T to degrees Kelvin 
and checks for out of range. If KR=1, T is checked for 
out of saturation range. If T is out of range, the pro- 
gram writes a note and continues. 

This routine is used frequently to insure that the proper regions and 
units are used in the calculations. It prints OUT OF RANGE messages to 
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warn the user that the results are extrapolated, and may be incorrect. 

It is the responsibility of the user to check these results for validity. 

FUNCTION SOLVE (XI, F,DF) 

This routine is a Newton-Raphson iteration given an initial estimate 
XI, the function F, and its derivative DF. If the convergence has 
not been attained in 100 iterations, the result and a note are 
written out. Calculation continues with the last value. 

SUBROUTINE ROOT (X0 ,X2 ,F0FX,FUNC ,X1) 

This routine solves for XI such that FUNC(X1)=F0FX, where XI 
lies between XO and X2. It includes a modified half -interval 
search technique and permits only one root between XO and X2. 

If a root has not been found in 100 iterations , the result and a 
note are written out , and calculation continues . 

SUBROUTINE R00TX - Same as SUBROUTINE ROOT 

This routine must be included to allow a function argument of R00TX 
to call ROOT. 

SUBROUTINE SPLINA(X ,Y ,NX,T ,NT ,YINT ,KFD,KERROR) 

This is a spline curve fit routine used for interpolation. 

SUBROUTINE P0LY(X,C0EF) 

This routine evaluates the polynomial at X, where f(x) is described 
by the coefficient array COEF. 

SUBROUTINE DENS (KU,T,P ,D,DL,DV,KR) 

This routine computes the density D, given the temperature T and 
pressure P, If KR=1, the saturated liquid and vapor densities, 

DL and DV, respectively, are computed as a function of T only. 

If either T or P is initialized to zero, the saturated densities 
are computed as functions of P or T, respectively, and the satu- 
ration temperature or pressure is returned in the initialized store. 

SUBROUTINE PRESS (KU,T,D,P,KR) 

This routine computes the pressure P, given temperature T and 
density D. If KR-1, P will be the saturation pressure computed 
as a function of T only by SUBROUTINE PSSS. While the equation 
of state is found in SUBROUTINE PRESS, it is also used in subsequent 
routines. The state equation used by GASP is: 
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P = RpT + A(T) p 2 + B(T) p 3 + C(T)p 4 + D(T)p 5 + n^p 6 


3 2 "" n 20^ ^2 

+ p [E(T) + p F(T)]e 


where if the fluid is helium, T 2 = T, otherwise T 2 = 1. 


n 3 n 4 n 5 n 21 
A(T) =n 1 T + n 2+ - + - I + - J + - r 


B(T) = n 22 T 


2 


+ n,T + n, 
6 7 




C(T) = n g T + n 10 


D(T) 


= n ll T + n 12 



E(T) 



1 


T 


2 


F(T) 




1 


T 


2 


The coefficients n-^ to are from Bender's original equation while 

n 2 i to n 24 have been added for carbon-monoxide, neon, and helium. 

SUBROUTINE TEMP (KU , P , D ,T ,KR) 


This routine computes the temperature T, given pressure P and 
density D. If KR=1, T will be the saturation value computed as 
a function of P only, by function TSS. 


FUNCTION DSF(DS) (A multiple entry routine) 


This is a function used with DENS to solve for density DS given 
temperature TS and pressure PS. 

DSF = P equation of state - P given 

ENTRY DDSF(DS) 

This routine calculates the derivative of the above func- 
tion and is used with DENS to solve for density DS given 
temperature TS and pressure PS 
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DDSF = 


8 (DSF) 

8p. 


SUBROUTINE PSSS(PSS) 

This routine computes the saturation pressure PSS given tempera- 
ture TS , where 

lo g 10 P - 3 ± + ^ + J3T.+ j 4 T 2 + J 5 T 3 + j 6 T 4 + j y T 5 

and the values to are found in SETUP. 

FUNCTION TSS(PS) 

This is a function used to compute saturation temperature given the 
pressure PS. 

FUNCTION TSSF(TSS) (A multiple entry routine) 

This function is used with FUNCTION TSS to solve for saturation tem- 
perature TSS given pressure PS. 

TSSF = Vapor pressure equation - log^CPS) 

ENTRY DTSSF (TSS) 

This routine provides the derivative of the above function 
and is used with function TSS to solve for saturation 
temperature TSS given the pressure PS 

DTSSF = 3 - ^" S - F - ~ 
a I 

FUNCTION TSF(TS) (A multiple entry routine) 

This function is used with TEMP to solve for the temperature TS, 
given pressure PS and density DS. 

TSF » Equation of state - PS 

ENTRY DTSF(TS) 

This routine calculates the derivative of the above func- 
tion and is used with TEMP to solve for temperature TS 
given pressure PS and density DS . 


DTSF = 
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SUBROUTINE TEMPPH (KU,P,H,T,D,DL,DV,KR) 

For a given pressure P and enthalpy H, this routine will return 
values of temperature T and density D. If KR=1, saturated 
liquid and vapor properties will be returned. 

SUBROUT INE TEMPP S(KU,P,S,T,D,DL,DV,KR) 

For a given pressure P and entropy S, this routine will return 
values of temperature T, and density D. If KR=1, saturated 
liquid and vapor properties will be returned. 

FUNCTION TSHF(TS) (A multiple entry routine) 

This is a function used in conjunction with TEMPPH. It 
obtains a trial value of DS using the given PS and a 
trial TS. Then it obtains a trial H which is compared 
to the input enthalpy within an iteration in TEMPPH. 

ENTRY TPSF(TS) 

This is a function used in conjunction with TEMPPS. It 
obtains a trial value of DS from the given PS and a 
trial TS. Then it finds a trial S which is compared 
to the input entropy within an iteration in TEMPPS. 

SUBROUTINE ENTH (KU ,KR,T,P ,D ,H,HL ,HV) 

This routine computes enthalpy H given the temperature T, pressure 
P, and density D. If KR=1, the saturated liquid and vapor enthal- 
pies HL and HV, respectively, are computed as a function of T only. 

SUBROUT INE ENT (KU , KR , T , P , D , S , SL , SV) 

This routine computes entropy S given temperature T, pressure P, 
and density D. If KR=1, the saturated liquid and vapor entropies 
SL and SV, respectively, are computed as a function of T only. 

FUNCTION HDINT (DS ,DSL) (A multiple entry routine) 

This routine computes the integral used in the enthalpy 
computation from density DSL to density DS. 
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ENTRY SDINT (DS ,DSL) 

This routine computes the entropy integral from the 
density DSL to the density DS 



FUNCTION HDINTF(DS) (A multiple entry routine) 

This routine evaluates the integrand of function HDINT 
where DS is the variable of integration. 



ENTRY SDINTF(DS) 

This routine evaluates the integrand of function SDINT 
where DS is the variable of integration. • 



FUNCTION HSS (PS ,DS) (A multiple entry routine) 

This function computes the enthalpy in the region KR=3 
(see table IV) or saturated vapor enthalpy (KR=1) , given 
pressure PS, density DS, and temperature TS . 
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and the values m^ to m^ are found in SETUP. 

ENTRY SSS (PS ,DS) 

This routine computes the entropy in region KR=3 or the 
saturated vapor entropy (KR=1) given the pressure PS , the 
density DS, and temperature TS. 
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SUBROUTINE HSLV(PS) (A multiple entry routine) 

This routine is used in conjunction with HSS to compute 
the saturated liquid and vapor enthalpies given pressure 
PS and temperature TS. HV comes from HSS, and 

HL = HV - T*HSSLVF 


ENTRY SSLV(PS) 

This routine is used with SSS to compute the saturated 
liquid and vapor entropies given pressure PS and tem- 
perature TS. SV comes from SSS, and 

SL = SV - HSSLVF 


FUNCTION HSSLVF ( PS) 

This is a function used with HSLV to compute saturated liquid en- 
thalpy from saturated vapor enthalpy or saturated liquid entropy 
from saturated vapor entropy given pressure PS and temperature TS. 

HSSLVF = ££ AV 


where 


P 

SUBROUTINE CPPRL (PS ,DS ,T,CPPART ,CVPART,KU,KR,KCP) 


This routine computes the specific heats at constant volume (CV) 
and constant pressure (CP). IF KR=1, the saturated liquid and 
vapor specific heats at constant volume, CVL or CW, and constant 
pressure CPL or CPV, respectively, are computed as requested by 
GASP. 


The specific heats are computed by 


CV - Cp Q - R - 



dp 
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T (3P / 3T) ^ 

Cp - CV + ~2 Op/ 3p) ^ 

provided the input density is less than p . , (p = critical 

density). Also see figure 36. “’' 1Ca *" 


Fluid 

P switch^c 

Helium 

2.4 

Methane 

2.5 

Neon 

1.0 

Nitrogen 

2.2 

Carbon monoxide 

2.4 

Oxygen 

2.4 

Argon 

2.4 

Carbon dioxide 

2.4 


If density is greater than these values, the specific heats are 
calculated using SUBROUTINE SPCHV. 

In addition, if T < T C ritical> th® 11 the specific heats for neon, 
carbon-monoxide, and helium will be computed by SPCHV. 

The partial derivatives (3P/3T)p and (3P/3p)T are available to 
the user in COMMON/ PARTLS/PTV,PDT. 

FUNCTION CPPRLF(DS) 

This routine evaluates the integral used in computing specific heats. 



SUBROUTINE PTRHO (D ,T) 


This routine evaluates two partial derivatives . The answers are 
returned in COMMON/PARTLS/PTV ,PDT . 


P 


PTV - 


3P 

3T 


PDT = 





14 


SUBROUTINE SPCHV(KU,KR,T,P,D,CV, CVL , CW) 


This routine computes the specific heat at constant volume CV 
given temperature T, pressure P, and density D. If KR=1, the 
saturated liquid or vapor specific heat, CVL or CVV, respectively, 
is computed as requested by GASP. 



See SUBROUTINE CPPRL. 


SUBROUTINE CVPS (KVP ,KR,CVS) 

Given temperature, pressure, and density, this routine is used to 
determine five values of internal energy U for a spline fit used 
in SPCHV to compute the specific heat at constant volume. 

SUBROUTINE SONIC (KU, KR, T ,D, GAMMA, C) 

This routine computes sonic velocity C given temperature T, 
density D, and the specific heat ratio GAMMA=CP/CV. 

SUBROUTINE VISC(KU,KR,T,D,MU) 


This routine computes the viscosity MU, given temperature T and 
density D. The dilute gas viscosity js computed as a simplified 
form of the Chapman-Enskog model by p = 0.2669 E - Av^aT/aa^fiy , 
where £2^ is the Leonard-Jones potential viscosity collision inte- 
gral, reference 16. The excess viscosity correlation is from Josci, 
Stiel, and Thodos , reference 13. 


* -4 1/4 

[(y - y )?xl0 ] = 0.1023 + 0.023364 p_ + 0.058533 p 


- 0.040758 + 0.0093324 p£ 

K K 


where y - y data has been verified in the references listed in 
table I. Here p R = p/p c and g = Tl/6y^-/2p2/3 

SUBROUTINE THERM (KU , KR , T ,P ,D,EXCESK,FK) 


This routine computes the thermal conductivity FK, given the tem- 
perature T and density D. FK represents the thermal conduc- 
tivity of the dilute gas plus the excess conductivity (X - X*) which 
is a function of density. The computation of the thermal conduc- 
tivity anomaly EXCESK (X R ) is based on the theory of Sengers and 
Keys, reference 43, which agrees with the reacting gas theory of 
Brokaw, reference 41. The techniques are compared and discussed in 
reference 44. 
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EXGESK will be computed for the density regime 0.3 < p/p c < 1.8. 
To obtain the thermal conductivity of a near-critical fluid, the 
user must add FK and EXCESK . 

The dilute gas conductivity and the excess conductivity are curve 
fits from the data in the references as listed in the fluid refer- 
ence table I* The Stiel and Thodos correlation of reference 36 is 
used where it is practical. Special procedures are followed for 
helium as indicated in the references of table I. 

SUBROUTINE SURF (KU , KR , T , SIGMA) 

Compute the surface tension SIGMA, given the temperature T, by a 
corresponding-state correlation suggested by Brock and Bird, ref- 
erence 39. # 

SUBROUTINE DGUES (TS ,TCR,DST) 

This routine provides a near critical density estimate used by 
subroutine DENS, 

FUNCTION CONZ (TEMP) 

This routine computes the thermal conductivity of helium for the 
dilute gas as a function of temperature by the technique of 
Roder , reference 55 . 

In closing, now that the reader is familiar with GASP, turn to 
page 3 1 , excise tables IV and V and discard the remainder of the paper. 
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APPENDIX B 

A LIST OF SUBROUTINE GASP 


$[ BFTC GASP Y DECK 

C 

C NASA LEWIS RESEARCH 

C VERSION 4/19/71 

c FOR INFORMATION —SEE A. BARON* R.C.HENDRI CKS* I.PELL ER 

C 

SUBROUT INE GASP! KS,KP*T,P *0*H*KR > 

C BENDER— S EQUATION OF STATE FOR THE FIVE GASES — N2 *CH4 * AR* 02* C02 

C ALSO EQUATIONS OF STATE FOR CO. NE * ANO HELIUM... 

C PRYDZ EQUATION FOR FLUORINE 

C 

C COMPUTF THE STATE RELATIONS ANO THERMODYNAMIC AND 

C TRANSPORT PRQPERl TES OF SPECIFIED FLUID GIVEN TEMPERATURE T, 

C PRFSSURE J>* DENSITY D, OR ENTHALPY H. STATE RELATIONS ARE 

C SPEC IF I EO BY KS. THERMODYNAMIC AND TRANSPORT PROPERTIES 

C ARE SPECIFIED BY KP. IF KR IS RETURNED OR SPECIFIED AS 1* 

C PROPERTIES ARE COMPUTEO AT SATURATION. 


C 


c 

c 

c 


c 

c 

c 


c 

f. 

c 


c 

c 

c 


c 

c 

c 


DIMENSION KP;ii32). KPC2I32)* KPC3 ( 32 I . KPC4 ( 32 ) 

COMMON/ PROPTY/KU.OL* DV »HL *HV* S* SL, SV*CV *CVL *CVV *CP* CPL »CPV * GAMMA* 
IGAMM AL.GAMMAV.C.CL*CVP.MU.MULtMUV.K*KLf KVtSIGMAtEXCLt EXCV tEXCESK 
REAL MU.MUL*MUV*K*KL.KV 

COMMON /CHECKS/DCH1 *DCH2*PCHl , PCH2 *PCH3 ,TCH1* TCH2* TCH3 *OST* TST* H 
1SCH1 . HSCH2 

DATA KPC1 /2*3*6.7* 10,11*14,15,18*19*22 *23,26,27*30,31*34*35,38, 
139.42.43*46*47, 50,51. 54*55*58,59.62,63/ 

DATA It PC? /4 *5*6*7* 12 , 13, 14, 15 , 20 *21 ,22 *23, 28 * 29 , 30 , 31 * 36 * 37* 38, 
139. 44, 45, 46* 47 *52, 53, 54*55*60*61. 62, 63/ 

DATA KPf-3 /8* 9* 10. 11 , 12 ♦ 13 , 14 , 15 * 24* 25* 26 . 2 7, 28 ,29 , 30 * 31 * 40 ,41 ,42 , 
143,44.4 5* 46* 55* 56 *57*58, 59*60,61 *62*63/ 

DATA KPC4 /16, 17, 18, 19. 20, 21*22. 23. 24*25*26*27*28, 29*30, 31, 
148*49,50*51,52*53,54.55.56*57*58*59,60*61 .62,63/ 

GO TO C 10. 20*30, 40.451 *KS 


COMPUTE OENSITY 


10 IF 1 KR • NF • l.AND.IP.LT.PCHl.OR. T.LT • TCHl WRITE (6,2511 
IF I KR.EQ. 1. AND. (P.LT.PCHl. AND. T.LT.TCH1) ) WRITEI6.252) 
IF ( P.GT.PCH3.QR .T.GT.TCH3) WRITEI6.253) 

CALL OENSI KU • T.P *D*DL *DV * KR ) 

GO TO 50 


COMPUTE PRESSURE 

70 IF (KR.NE.1.AN0.IT.LT.TCH1.QR.D.LT.DCHU > WRITE(6,251) 
IF < KR.EQ. 1. ANO. (T.LT. TCHl. ANO. D.LT.DCHin WRITE<6,252> 
IF ( T.GT.TCH3.0R.D.GT.0CH2I WRITEI6.253) 

CALL PR ESS<KU,T*D*P,KR I 
GO TO 50 

COMPUTE TEMPERATURE 

30 IF (KR.EQ. 1. AND. (P.LT.PCHl. ANO. O.LT.OCHDI WRITE(6,251) 
IF ( KR.NE. 1. AND. (P.LT.PCHl. OR. D. LT. DCH1 1) WRITE (6 *2521 
IF (P.GT.PCW3.QR.D.GT.DCH2J WRITE<6,253) 

CALL TFMP ( KU* P*0*T.KRi 
GO TO 50 


COMPUIfc TEMPERATURE ANO OENSITY GIVEN PRESSURE AND ENTHALPY 

40 IF ( P.LT.PCHl. OR. P.GT.PCH3) WRITE(6,254I 
CALL TEMP PHI KU.P.H.T.OfDL.DV.KRi 
GO TO 50 


COMPUTF TEMPERATURE ANO DENSITY GIVEN PRESSURE AND ENTROPY 

45 IF ( P .LT. P CHI .OR. P .GT. PCH3 ) WRITE (6,254 I 
CALL TEMPPS ( KU .P . S , T , D , DL *□ V , KR ) 

50 IF (M0DIKP.2)) 60*70, 6C 
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G 

C COMPUTE ENTHALPY 

C 

60 CALL EN THI KU, KR ,T* P ,D,H«HL«HVl 
70 00 80 1=1*32 

IF ( KP-KPC 1 ( I ) ) 1 10 1 100 * 80 
80 COMTINUF 
G U TO 110 
C 

C COMPUTE ENTROPY 

C 

100 CALL ENT!KU,KR,T,P,D,S,SL,SV) 

110 00 120 1=1*32 

IFUP-KPC2(IU 140*130* 120 
120 CONTINUE 
GO TO 140 
130 KCP=0 
C 

C COMPUTE SPECIFIC HEATS ANO SONIC VELOCITY 

C 

IF { KR « NE • 1 ) GO TO 20C 
CALL UENS<KU,T,P,D,DL,DV,1) 

CALL CP PRL ( P* DL * T • CP L»CVL ,KU, KR * KCP I 
CALL CP PRL ( P* OV « T » CP V *CVV *KU , KR , KC P) 

GAMMAL=CPL/CVL 
GAMM A V = CP V/CVV 

CA LL SONI C ( K J • KR , T , OL , GAM MAL . CL ) 

CALL SONIC !KU,KR,T,DV, GAMMAV* CVP I 
GO TO 140 

200 CALL CPPRL(P.D*T,CP*CV*KU*KR*KCP) 

GAMM A=CP/C V 

CALL SONI C ( KU* KR , T * D*GAMMA*C ) 

140 00 150 1=1,32 

IF IKP-KPC3UII 170,160*150 
150 CONTINUE 
GO TO 170 
C 

C COMPUTE VISCOSITY 

C 

160 IF ( KR • NF • 1 ) GO TO 165 

CALL OENS { KU, T ,P *D, OL *DV, 1 ! 

CALL VISCI KU,Kf**T,DL,MUD 
CALL V I SCI KU*KR,T*DV*MUV) 

GO TO 170 

165 CALL VCSC(KU*KR*T,0,MU) 

1 70 on 175 1=1*32 

IFUP-KPC4UM 190,180*175 
175 CONTINUE 
GU TO 190 
C 

C COMPUTE THERMAL CONDUCTIVITY 

C NOTE— FRTZEN VALUE AVAILABLE IN K, KL * K Y 

C — REACTING CONDUCTIVITY RETURNED IN EXCESK, EXCL.EXCV 

C 

1 HO IF ( KR » NE • 1 1 GO TO 220 

CALL OENS I KUtT*P*D*OL,DVtl» 

CALL THERM ! KU, KR ,P * T ,DL • EXCL * KL J 
CALL THERM (<U , KR, P , T ,DV * EXC V * KV) 

GO TfJ 190 

220 CALL THERM < KU* KR. P ♦ T ,0 , EXCE SK * K I 
190 IFKP-32) 230,240*240 
C 

C COMPUTE SURFACE TENSION 

C 

240 CALI SURF ( KU, KR • T, S IGMA I 
230 RETURN 

261 FORMAT! 1H0.45H OUT Of. RANGE -LOW- SEE ERROR L I ST-CALC. CUNT 

252 FORMAT! 1H0,45H OUT OF RANGE- SAT- SEE ERROR L I S T-CALC • C ONT 

253 FORMAT! 1H0.45H OUT OF RANGE-Hl -SEE ERROR L I S T-CALC. C ONT 

264 FORMAT! 1H0.45H OUT OF RANGE- P -SEE ERROR L I S T-CALC. CONT 

FNI) 
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APPENDIX C 
SYMBOL LIST 


Fortran 

symbol 


c 


Cp 

cp 0 


Cv 


P 


P 


t 


1 


R 


G 

CL 

CVP 

CP 

CPOCOR 

CPL 

CPV 

CV 

CVL 

CVV 

D 

DL 

DS 

DSL 


DV 

EXCL 

EXCV 

EXCESK 


sonic velocity, cm/sec 

sonic velocity, saturated liquid, cm/sec 
sonic velocity, saturated gas, cm/ sec 
specific heat at constant pressure, J/ (g) (K) 
specific heat at "zero' 1 pressure, J/(g)(K) 
specific heat units correction factor 
saturated liquid specific heat, Cp, J/(g)(K) 
saturated vapor specific heat, Cp, J/(g)(K) 
specific heat at constant volume, J/(g)(K) 
saturated liquid specific heat, Cv, J/(g)(K) 
saturated vapor specific heat, Cv, J/(g)(K) 
density, g/cm^ 

3 

density, saturated liquid, g/cm 

3 

density used internal to the program, gm/cm 

density, saturated liquid, used internal to the 
program, gm/cm^ 

3 

triple point density, g/cm 

3 

density, saturated vapor, g/cm 

reacting conductivity, saturated liquid, w/cm-k 

reacting conductivity, saturated vapor, w/cm-k 

i 

reacting conductivity, w/cm-K 


Symbols used in the individual subroutines are identified in the 
work statement of that subroutine (see appendix A) . 
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Y GAMMA 

GAMMAL 
GAMMAV 
H H 


H 


0 


3 


1-7 


k 


* 


k 


HV 

HL 

set 


K 

KCP 


KL 

KP 

KR 

KS 

KU 

KV 

7 7L 

m l-5 

y MU 

MUL 
MUV 

n l-24 

P P 


ratio of specific heats, Cp/Cv 

ratio of specific heats , saturated liquid 

ratio of specific heats, saturated vapor 

enthalpy, J/g 

reference enthalpy, J/g 

enthalpy, saturated vapor, J/g 

enthalpy, saturated liquid, J/g 

set of constants, equation of vapor pressure curve 
thermal conductivity (p -* 0) W/ (cm) (K) 
thermal conductivity, W/ (cm) (K) 

region delimiter used in the numerical calculation 
of CV 

thermal conductivity, saturated liquid, W/ (cm) (K) 
thermodynamic and transport properties specification 
thermodynamic region specification 
state relation specification 
units specification 

thermal conductivity, saturated vapor, W/(cm)(K) 
molecular weight 

set of constants, equation of specific heat at "zero" 
pressure 

dynamic viscosity, saturated liquid, g/(cm) (sec) 
dynamic viscosity, saturated liquid, g/ (cm) (sec) 
dynamic viscosity, saturated vapor, g/ (cm) (sec) 
set of constants, equation of state 
pressure, atm 
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P 

c 

P 

max 

P . 
min 


^switch 

S 


a 

a 

S 


d 

0 


T 


T 


c 

T 

max 


PC 


PS 

PSS 

R 

RHOC 


S 

SIGMA 


SL 

SV 

T 


TC 


TS 

TSS 


t 

U - H 


P 

P 



P 

pRT 


pressure at the thermodynamic critical point 

2 

upper pressure limit, MN/m 

2 

lower pressure limit, MN/m 

2 

pressure used internal to the program, MN/m 

saturation pressure used internal to the program, 
MN/m2 

gas constant, J/(g)(K) 

3 

density at the thermodynamic critical point, gm/cm 

density where calculation of CV change from numeric 
to analytic 

entropy, J/(g)(K) 

surface tension, dyne/cm 

hard sphere collision diameter, angstroms 

reference entropy, J/(g)(K) 

entropy, saturated liquid, J/(g)(K) 

entropy, saturated vapor, J/(g) (K) 

temperature, K 

reference temperature, K 

temperature of the thermodynamic critical point, K 

upper temperature limit, K 

temperature used internal to the program, K 

saturation temperature computed by FUNCTION TSS, K 

triple point temperature, K 

internal energy, J/gm 

specific volume, cm /gm 

compressibility factor 

Leonard-Jones collision integral-viscosity 
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APPENDIX D 

COMPACT FLOW CHART FOR GASP 


CALL DENS 

(KU,T,P,D,DL,DV,KR) 


ENTER 
OUTINE GASP 


(KS,KPT,PD,H,KR)) 



CALL PRESS 
(KU,T,D,P,KR) 



U> 


CALL TEMFPH 
(KU,P,H,T,DL,DV,KR ) 



CALL TEMP 
(KU,P,D,T,KR) 
^ 




Tests are so arranged, but 
not indicated here, such 
that i.e., KP= 1+4=5 
will return enthalpy and 

<^KP=1 

N0| 

Jk 

specific heats (CP, CV77" "their 
ratio = CP/CV, and sonic 

<KP=2 

velocity. See, also, the 
operations sheet. 

noJ 


<^KIfc4 


NOj 


YES 

CALL ENTH 


(KU,KR,T,F,D,H,HL,HV) 


. YES 

CALL. ENT 


(KU,KR,T,P,D,S,SL,SV) 



YES 


CALL CPPRL 
(F,D,T,CP,CV,KU,KR,KLP) 


CALL DENS 
(KU,T,P,D,DL,DV,1) 


CALL CPPRL(P,DL,T,CPL,CVL,KU,KR,KCP) 
CALL CPPRL(P,DV,T,CPV,CVV,KU,KR,KCP) 


GAMMA- CP/CV 


|gammal= cpl/cvl j 
|gammav= CPV/CW I 

3: 


CALL SONIC (KU,KR,T,DL,GAMMAL, CL) 
CALL SONIC (KU,KR,T,DV, GAMMA V,C VP )| 


I 1 

CALL SONIC 
(KU, KR,T,D, GAMMA, C) 




CALL DENS 

^ yes ( 


_ YES 

(KU,T,P,D,DL,DV,KR) 


V KR 1 y 




i - 

Tno 


CALL VISC (KU,KR,T,D, MUL ) 
CALL VISC (KU,KR,T,D, MUV ) 


CALL VISC 
(KU,KR,T,D,MU) 


r- 


1 



RETURN \ 

(answers to user' s programy 
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r ENTER 

SUBROUT INE DENS (KU,T,P,D,DL,DV,KR ) , 


|ps-pcheck(ku,kr,p)H — *° 


TS=TSS (PS ) 
T-TS 


1 TS^TCHECK { KU , KR , T )\ 


[CALL PSSS(PS) 
P=PS 


INITIALIZE DS 


TS=TCHECK(KU,KR,T) 
PS= PCHECK(KU,KR,P) 


I *DSL a SOLVE ( DS , DSF , DDSF )1 


INITIALIZE DS 


:heck for\ 
REGION NUMBER 
i.e. set KR-2 
\or KR=3 ^ 


INITIALIZE EST 




♦See APPENDIX A for details of subroutine or function specified. 
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♦See APPENDIX A for details of subroutine or function specified. 


TEMP uses : 

PC HECK 
DC HECK 
TSS(PS) 
DENS 
SOLVE 



♦See APPENDIX A for details of subroutine or function specified. 
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*CALL HSLV(PS) to 
compute HSL, HSV 



NO 

r ■ ■■ 

DS*DCHECK ( KU , D)j 


TCHECK 
PCHECK 
DC HECK 
HSLV 
PSSS 


HL= HSL* HC ONV ( KU 
HV=HSV*HCONV(KU 


-1CALL PSSS(PSS)! 


CALL HSLF(PSS) to 
compute HSL, DSL 
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/ ENTER > 

( SUBROUTINE TEMPPH{KU,P,H,T,D,DL,DV,KR) 

^SUBROUTINE TEMP PS (KU , P , H , T , D , DL , DV, KR )/ 


♦TEMP PS is identical to 
TEMPPH except that H 
representing enthalpy is 
everywhere replaced by S 
to indicate entropy and 
ENTH is changed to ENT 
and TSHF to TPSF . 


CALL TEMP to get 
SATURATION TS 


CALL ENTH to get HSL, HSV 
j* CALL ENT to get SSL, SSV 


CALL DENS to 
get DSL,DSV 




|PS= PCHECK (KU,KR,P)| 




HS«H/HCONV(KU) 

ss**s/sconv(ku ) 



YES \ NO 

" 0>S<PCH2> ■ 


DETERMINE approximate REGION 
NUMBER on basis of: HSL, HSV, 

and PS 


DETERMINE approximate REGION 
NUMBER on basis of: SSL, SSV, 

and PS 


YES 




NO 


TEMPPH uses: 

PC HECK 

HCONV 

TEMP 

ENTH 

DENS 

ROOTX 

TSHF 


TEMPPS uses: 

PCHECK 

SCONV 

TEMP 

ENT 

DENS 


ROOTX 

TPSF 


►(INITIALIZE TS1, TS2| 



CALL ROOTX (TS1. 
♦CALL ROOTX (TSlj 

iTS2,HS,TSHF,TS) 
,TS2,SS,TPSF,TS ) 

\ 


[CALL DENS 

to get Ds( 


D*DS »DC ONV ( KU) | 


DL=DSL*DCONV(KU) 

DV»DSV*DCONV(KU) 


”T’ 

(return) 


VERIFY REGION NUMBER from 
given PS and calculated TS 
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VISC uses: 

TCHECK 

DCHECK 

POLY 




THERM uses: 

TCHECK 

DCHECK 

POLY 

CONZ 


(return) 
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SURF uses: 
TCHECK 



CPPRL uses: 

CPPRLF 

3PCHV 

PTRHO 




































































TABLE II. - COEFFICIENT MATRIX FOR SUBROUTINE SETUP (NAMGAS) 


29 


O 

o 


u 

< 


w 

o 


rii~i n n r- 

IM yJ- 

in 

If- li yj- in 4) f- N- CD fM cn 

^1 

Oo © o o 

O O 

© 

OOOOOOOOOtJ 

<J o 

+ + + + + 

+ + 

+ 

+ 4-U + + ++ + + + 

+ + 

LU LL' UJ LiJ U-i 

U! LJ 

LLl 

UiUi UJ UJ UJ U.: Ui U.1 U' ui 

UJ UI 

4- in PP) CM QD 

in ^1 

43 

4) in f\| CO P- 0 © © 43 rH 

43 00 

N 4 >J 4j 4 

03 CM 

43 

OOCVlN-43'ta343fMCMin 

CM <4- 

OJ CM 4) 4) in 

r-l in 

OD 

ytnONHNinttJN 

oo N- 

pm r IN N N 

© m 

43 

»+ 00 43 O r-l © AJ N* Off 

o in 

1— N- i la m 

m in 

CO 

njN<rrti4 , in'0tnNci 

© m 

N 4 4 03 ff- 

nj 4) 


fON oin O' N-4- coin cn 
4)00'SMfiCOH4)00 

N- CM 

C 4 rtN H 

O -4- 


r-l CM 

OJ CO CM y+ r-l 

CM r-. 

r-l 

N N N Ifl r-l rl yf rl 4' rt 

cm cn 

do o o o o 
t i 

oooood 

i 

©o ooooooooooo 

II II 

00 

1 1 

OONifyD 

O CM 

O 

Pi N © H CO O (VI rl O’ 

cm cn 

OO o o o 

OO 

O 

OOOrHOrHr-lOO 

oo 

+ + + + + 

+ + 

+ 

4 . + + + + + + + + 

+ + 

II 1 1 1 1 III II 1 III 

UJ UJ 

IU 

UiLUUJUJUJUJUJUJUJ 

UJ UJ 

N 4 rt r-i in 

C N- 

r- 

N-N“N-©ms|-cny4'© 

CO 00 

O yj- cn -4- O' 

rH © 

in 

4"4- N fO Mn (V -4- n 

£S 

cd in in © o 

n- -4- 

H 

NHNHC'H4-4rt 

CM P- rH cn © 

03 rH 

in 

(OHiONN 4)4)Nin 

CM CM 

Qy -yj" yj 1 rH PM 

O m 

43 

'0«0 4'N414-C0rtN 

N- CO 

4J in CM N- rH 

rO rH 

<n 

m O © >4- N“ in CM 43 N- 

rH vN 

Ct 4) O 4) rt 

in © 

cn 

rH 00 fp| yO 4)N4iC 

cn cm 


(\J m 00 H 


hoO^M 
o o o o c 
+ + + I I 
Uj u.i UJ lu ut 
rH cm N- c in 
O0'(M<ON 
00 ^ r-. «o 
H O' h o o 
© O' O O' O' 
in O' f" i-i h 

cn m nt © cm 

W O' M N H 

• • • a a 

0.0 o o o 

I I 

S cm yj- n- O 
O O O rH 
+ 11*1 
LU UJ UJ UJ UJ 

oo oo c* O' rvi 
co cm in m 
coow'tcn 
N- in o> co © 

OCOO"f N 
NfO sO «1->D 
® O cn rH ** 
»-i cm »-i in oo 



I I 


CM 4 
OO 
+ + 

UJ UJ 

cm rn 

rH in 
rH © 
cn CO 

-* rt 

0 r~ 

in so 

O' CM 

• .» 

0 o 

1 

rvi an < 


-h cm 4 N> O 
O o r» C- rH 
+ 1111 
UJ UJ UJ UJ UJ 
3 n- in « 


P 


in r-l s0 o 
m in o >o ■o 
n sD n m in 

ON<MftlM 

«tONHM 

HSHh-JI 

rlJlNmiNI 

a a a a a 

oooro 

I > 


o © 
o in 
n- n- 
<7* ao 
n- « 
© cm 

O' O' 
CM 4 


I 

fO © 
O O 
+ + 

S3 

in o 

00 rH 

4 4 
CM 4 
cn rH 


w — CM 4 aKflOlfllflfyO'HlSO rH CM 
OOOUOOOOOOOrHCrHrHO 
+ I + ++<* + 4-+ + + + + + + + 

UJUJUJlUUJUJUJlUUJUJUJUJUJUJUJUJ 

h*NOHaiOM)'<MpiMn(S)MO 

rt^aimooDH^ino'rnir'aiNoo 

'J , 0^aJN(C>tN«3in0r<fCNlJ'0 

rvcomNONtOfC'O^HMnjHiflci 

cciwoinMe^iPanft^r-iPrtNO 


i 


I 


J H in 00 O' N IT 00 ytallM^NaJ^lflal 


OCSQOQCOOOCOCOOCOCCCCOCCC 
til II 11)1 


O C 
+ + 

LL LL 

a m 

a r- 

rH CM 

in in 
cm n 
oc <n 
- 4” 
n rH 


cn in co 

$55 

U. UJ LL' 

-H rH J 5 

n- a* «n 

NfC + 

O' vC in 
\J) (M OC 
U) 4 PI 

-» o o 

O + H 


O' H (f U) rvj 4 IM ^in«famff'(M(CiN 

coco ©r» a cco^-ic^^c 

+ 4- + + ++ + ++ + + + + + + + 

UJ LU UJ LU IL LU UJ [XjLLUJ LULL'LU UJ UJ LL 

n sir ao'mocHino'^vONMno 
4 r-i r-C 0* Cn H (J> 4 MP C CO O' O 
cm 00 © CM cC N- © CinOnHinitriC 
<ooMfl'4inffli-.4 , inirh-cflo'Na 
m m a « 4 + 4 4ru O’ nun o> o> -oo 

in M CC J5 0C r* If H4’U)(T'r>Orl4C 

>0 oc n- oo rvewoc O'yOinn-cMO'ramp- 

PV HI vC CM rH © rH UlrtTrlr, 4 hM<I 


I 

© O' 

o o 

is 

© 45 

en CT 
4 in 
co cv 
O' O 

in O' 

cm 4 

d d odd oo 

ii I i 

(VIA CN4 4ia 

OC OOC CD 


rim On4 

oo o o o 

+ + + I I 

. 1 1 . 1 1 It I 1 1 1 uj 

© CO •— I O' O' 
CM 4) as CM in 
+ C 4 H-. 
ff> 4) N (M N 
©CO sn 00 
CO O O © «“i 
'CONN'D 
© rv cn © in 


-.ill 
lu Uj UJ Uj Uj 
o rv <fi O' a 
r-nv in cn m 
Nhh 4 O' 
ico'riea: 
n- 4 in ro 4 
Mf CC dm 
ONC N4 
N- N- 4 « 


j <? a n a a onocomConcc one 

II II 


Or-t<Mc04Lntot^aocnOi-iCvito-4 
rH CM tO 4 If 3 CO t - -- GO CJJ i — ] c — t i — I H H i — 1 1 — i i — I i — t i — l CM CVl CM CVI OJ 

Pf; fl C! C C C C C C fl c c c c c c c c c c c c c c 


lL uj 
© m 

O' rt 
O' O' 
r** 43 

o* n- 

o O' 

CM O' 

cn cm . 


i c r> c c n 


oo o o o 


"i in i 


OO o 

+ I I 

LU Ut UJ 

cm c j cn 

C'fM CM 

© cm in 

CMm CM 
O' m cn 
(Cm u 
C'O © 
<-* in cm 


oo © o © o 
1 I 

ooooooooooooooooooo 
1 III 

O O o © © © © 
It 1 1 

or- 

i 

OON4’lfl 

rH (M 4 O N 

CM 43 03 

N- CD 

CM 

cm cn m 

rH cn © cm m in oo 

rH 

o © o c> o 

+ + + + + 

O O C» O C' 

| + 4* + + 

tlt’O 

+ + + 

5$ 

¥ 

T S* 

O e. p c OOC' 

+ + + 1 1 1 1 

© 

+ ■ 

UJ Ui UJ LLl LU 

UJ LU UJ UJ UJ 

UJ UJ UJ 

LU UJ 

LU 

Ui UJ UJ 

UJ UJ LU lu UJ UJ UJ 

UJ 

in N- m co ro 

C' 43 in yQ IM 

© © 4) 

CO rH 

cn 


in N- m © pm -4- m 

o 

rtN rt CO oo 

co r- cm n o 

00 CM © 

© O 

m 

cn cm m 

© cn 4- 1 ^- cn »h in 

o 

•4* cm cn cn in 

N- cn n- cm oo 

NNrt 

O N“ 

H 

"4" © *n 

O'® ©O' CO HO 

© 

in © n- co in 

NN 4 Hin 

4 , '0 4 

03 m 

© 

cn co oo 

cm © © © m © cn 

© 

00 rH © C3 

N- m O © 43 

© m y* 

CM © 

© 

mom 

CM© -4- CD CM CM N- 

© 

H00 vO n N 

pn co m © in 

© CD & 

© © 

in 

m co © 

4'N © © © CM N- 

CM 

rH CM rH in rH 

CM —4 tn 4) 43 

rH fCl O 

© © 

cn 

CM © 4” 

N *H © © 4* N* cn 

© 

4- -4- -4- rH CM 

in cn *h co ® 

yt CM rH 

rH N- 

rH 

HN# 

H H 4 H H © 

rH 


a a a a a 

a a a a a 

a a a a 

a a 

a a a 

a a • a • a a 

a a 

0 

0 

0* 

0- 

0 

9 

o o © OO 

O OD © O 

0 

*0 

0 

0 

Cl © 

© O © 

OD O O C> O O 

O© 


HaH in OO 
On OfJ 

♦ III 

UJ UJ UJ UJ 

O o «-iO 

oo oo 
OO OO 
04- no 

CM pyJ r-l m 

in cm oo rH 
Hin y0 cn 
lAH ON 

* • • • 4 

oo o o o 


fci 

o © 

o 

© 

o © 

O © 

O 

<■5 O 

o 

C5C3 

r3 o 

pH 

rt 

rH rH © 


55 
- o 

Cl o 

Q rt O 1*3 rH 

E 

o 

r> 

+ + 

+ 

+ 

4* + 

) + 

+ 

+ + 

+ 

+ + 

+ + 

+ 

+ 

+ + + 



+ 1111 


+ 


UJ UJ 

UJ 

UJ 

© © 

© © 

© 

© © 

© 

©© 

© © 

© 

© 

© © © 


© © 

UJ UJ © © © 

£ 

© 

a 

© o 

© 

o 

4 © 

rH rt 

CO 

ON- 

4 

rH PM 

00 t- 

CM © 

©CO © 



00 rH 

© CM O' N- C(3 

< 

rt 

O CM 

H 

N- 

4 © 

N* CM 

4 

PM CM 

in 

4 rH 

in N- 

in m 

© © © 



4 cn 

m © 4 ® cm 


O 


© © 

pH 

4 

m © 

© 4 

CM 

in cm 


pH rH 

© m 

N- 

rH 

co in o 



4 cn 

O CM in N- CM 

Ej 

o 

Eh 

cm m 

m 

N- 

© © 

CM rH 

CO 

N- m 

CM 

4 CM 

CM ® 

O oo 

CM © O 



4 © 

NQNrtN 

< 

© 


rH CM 

pn 

N- 

o © 

© © 

4 

© 4 

N* 

4 4 

m m 

4 m 

rt NO 


■ B 

I. rH © 

N. in in cm © 

K 

o 


© © 

N- 

N» 

4 © 

CM rH 

4 

t— © 

in 

in rH 

in in 

m 

OO 

©40 


I 1 

oo m 

Q cn © rH © 


a 

& 

O © 

H 

N» 

CM rl 

4 CM 

pH 

N- © 

pH 

© 00 

CM rH 

co m 

m r- in 


w 

cn n- 

CO oo o © e> 


m 

W 

CM rH 

© 

rH 

© cn 

4 © 

rH 

4 rH 

CM 

rH CM 

© © 

rH in 

rH © cn 


rH LA 

Hrl H rt ITl 

o 

CM 

• a 



a a 





a a 


a 

a 

, a a a • 

O OO © 
1 



i a a 

oo 

j 1 


H 

• a 


0 

0 

0 

1 

o 

0 O 

1 

OO 

1 

O 

uo 

1 

© 

1 

0 

-0 

o o 
1 

0 

1 

© 

1 

0 

0 

0 

, 

o© o OO 
1 1 1 

fxi 
s H 

OO 


o© 

© 

4 

N- © 

CM CM 

in 

o cn 

© 

CM CM 

N- © 

rH 

N- 

© rH rH 



1 rH IH5 

© PM 4 © © 

1 o 
w 

cm m 


© <-» 

o 

O 

© O 

O O 

o 

oo 

o 

o© 

o © 

rH 

O 

rH rH rt 



©o 

o © o © o 

pH 

oo 


+ + 

+ 

+ 

+ + 

1 + 

+ 

+ + 

+ 

+ + 

+ + 

+ 

+ 

H* + + 



. + + 

+ 1111 

CO 

+ 1 


UJ UJ 

UJ 

UJ 

© © 

© © 

© 

© © 

© 

© © 

© © 

© 

© 

© © © 



- © © 

© LU © © © 


UJ © 


© N» 

© 

4 

in 4* 

© CM 

© 

o © 

o 

rH © 

o o 

m 

© 

© N- O 



00 CM 

N- rH in 4 in 


O'© 


© N- 

r» 

4 

4 © 

© rH 

Q 

CD 4 

o 

N- N- 

eo 4 

© 

in 

CD © O 



rH CM 

m ® m « in 


<a © 


-H O 

© 

N 

CO rH 

cjy rH 

CM 

in 4 

CO 

N* O 

m cm 

CM 

© 

IM rH O 



4 © 

N. © © <n rH 


rH © 


CO rH 

© 


rH 4 

4 CM 

CM 

rH CM 

N- 

© nt 

rH 00 

© 

© 

4 N- O 



4 © 

O cm © 4 cn 


m o 


© rH 

N- 

© 

© N» 

cn N* 

© 

CM 4 

cn 

rH fPI 

cn © 

cn 

00 

4NO 



1 O CM 

© PM rH © CVJ 


4 rH 


© © 

O 

© 

O O 

rH rH 

© 

N- CM 

rH 

4© 

CM rH 

© 

© 

© © O 



in© 

CM CM m o 4 


'H O' 


in 4 

4 

in 

o © 

o m 
4 © 

O 

© © 

© 

© © 

© © 

© 

© 

cn rH 4 



frH in 

© © © o in 


0»© 


nj m 

rH 

CM 

rH rH 

rH 

© CM 

rH 

CM N- 

rH 4 

cn 

in 

rH © in 



in cm 

PM CO CM rH rH 


cm in 


a • 

a 

a 

a a 

a a 

a 

• a 

a 

a a 

• a 

a 

a 

1 « « • 

• • • 


a a 



a a 


oo 

O 

O 

O O 

OO 

O 

a- 

0 

o 

0 

0 

O © 

O 

© 

0 

0 

0 

0 

0 

0 

0 


O O 

0 

0« 

0 

0 

0 


OO 


r I r> r>, r-l 

oo o o 

+ + + + 

UI UI UJ U‘ 

oc o a 

OO <J o 
© 45 O O 
4-tMno 
4>n P's O 

Ntfi OO 

4 CM m rl 

od d d 


NrlfOiH 

coco 
+ + + + 
UJ LU UJ UJ 
o m o O 

04) on 

OP- a o 
Offl N- O 
4N CO O 

w~i O' cn o 
oo r- m o 

4)4 flf* 

a • • • 

oo oo 

NH Nat 

op O O 
UJ LU UJ UJ 

-Oh o © 
00 4 0 
ON NO 

ocion 
O' 4 4 © 
CJirtrto 
N4 NO 
CM CM O' *H 
a a a a 

oo o o 


oo 

+ + 
UJ UJ 
rt O' 


n N o 

CM rH 

pn rH 

O a rH 

OO 

o o 

1 1 1 

+ + 

+ 1 

© © © 

©© 

© © 

« CM © 

ON 

a co 

cn CM 

a © 

m ac 

o m O' 

OH 

00 CD 

4 © CM 

4 rH 

go cn 

O' c o 

© in 

rH 0* 

cn O rH 

cncn 

rH © 

CM © 4 

n- cn 

cn in 

O' CM rH 

N- CM 

cm m 

• • a 

a a 

a a 

COO 
1 1 

O© 

O C 

4 © O' 

-O rH 

n o 

O O O 

o o o u 


+ I 

Uj UJ 

mo 




+ 1111 

UJ XI UJ U. UJ 
mm rt N >. 
CM CM o a o 
4)0 OC -0 o 
tlfl n N m 
ON CO 4 O' 
© cv N- 4 m 
CN o >C rr t\| 
n- cn in <m cm 

a a # a * 

or: n n n 
I I 


i-l OJ tO LO 

X E S e E* 


c n 
O M- 
G y + 

CM 45 

m m 

rH N- 

e pm 

HNrlrt 

a • a a 


pn r-l 

o o 

+ I 
UJ UJ 
N- O 
r-« O 
i M- O 
i in r- 

l 45 45 
I 4) O 

i cm m 


a a • a 


O O 

o o 

CM CM 

m in 
O' aM 
r- as 


oo < 
* + 

UJ UJ 

oo 

or*- 

on 

'en 

GON- 

cm in 
N-N* 
oo CM 

OO 


0-4-0 0 
OOtON 
0 41 On 

§ CM CM a-l 

44IP 
O PM 4) CM 

o f— • in ^ 

O' m- cn m 

a a a a 

OO o o 


O in o 
O m O J s 
O NN4 

N rl O rl 

o m “M T> 
rtir, mo 

rH -n o o 

r-l cr> cu cvj 

a t a a 

0 O U U 

g 

o 

1 o 

O O Ocu 
f*r/)KO 


CM O 
cn o 

rH O 


rlm-IO 

n o G n ^ 

+ + I + sO 

LU OJ UJ UJ o 

O O O o o 

O O O Q 43 

O O O O I 

O o e» o c/3 

C O O C3 o 

o n rt 3 

N fi O' rl 

in 4) 4> PM 

a a a a 
0 0 3 0 


0 OO 

1 + + 

UJ UJ UJ 
ODC 

c> o c - 

uico 
CM 4) O 
m CP N- 
+H O' CM 
O yrf - O 
rH O IM 
. • • m 

ono 


cm cn pn 
O O O 
+ + + 
LU UJ UJ 
C'C. o 

o c< o 

CJOO 
O rH a 
M* O' o 
HN o 
ao m o 

4 )H 4 ) 

a a .a 

o r. c- 


o o 

4* + 
UJ UJ 

o c 
o c- 
oo 

Oft 

i^ r- 
O* M- 
O' CM 
CM r-l 


V? 

UJ UJ 

oo 
a e> 
in o 

CM N~ 
cn m 
—i in 

O 45 
WIM 


CM cm nj 
UJ uj LU 

o o o 
oo& 
© oo 
© © © 
O' <4“ a 
CM in M- 
O -+ -+ 
N N 4 


|f| OO 

UJ UJ UJ 
o o o 
o o © 
o oo 
o oc* 
o o © 
O cn 45 
© cn pn 
4) <# O 


oooooo OO 


r-laH CM 

on O 

I + + 

LU UJ UJ 

oo o 
o o o 
oo o 
cm in a* 

rn pr> 45 
H OD 41 
OrtO 
hn in 

a a a 

OO o 


oom 4 OH 
DOG DO 
+■ + + + + 
UJ UJ UJ UJ UJ 

o o o oo 
o o o oo 
O O © oo 

41 H ri OP) 

INNOOrt 

45-4t-©<NN- 

r-l O O 4) rH 

NrtH 4 h 

• a a a a 

OO O OO 


HH CM 

on o 

I + + 

Ut UJ UJ 

o© o 

OO O 
in© o 
cm© o 
in m 45 
*H 43 41 

O 00 o 
in 

a a a 


CNIrt 4 
O © © 
+ + + 
UJ UJ UJ 

©on 
OO o 
oo o 
o o o 

CO N- O 
N- O O 

m m o 

CO rH H 


O © 

4 + 

UJ UJ 

oo 

OO 
o o 
o« 


oo oooooo 


o© 

OO 

m© 

CM© 

m m 

-H 00 

OO 


CM CM cn m 
OOC© 
+ + + + 

UJ UJ UJ UJ 

o oo o 
o oo o 

O N- O O 

O O 00 o 

O tft N O 
Ort4C 
4‘4 , IOO 
m in —I m 

a a a a 

o o o o 


LU UJ 

o © 

O O 

o o 
o © 
m o 

CM rH 

in cn 

•4 rH 

a a 

a o 


00 i 

1 + 

+J uj i 
OQ 1 
a a t 
oo * 

CM© l 
IrtN I 


I o o o 
i o o o 
I o o o 
too© 
- N- n <n 

'oSSSS 

. rH ft) in 4 ) r-l 


UJ LU UJ 

o o o 
oca 
o o o 
© o a 
o in n 

DON 
O rH rH 
Cn rH 


ccddoddd 


U- LU 

e C 
o o 
a © 
© © 
N— N- 
45 CV 
rH 45 

rH 4- 


cv cm cn <n 
o n o o 
+ + + + 

LL LL 1 LU UJ 

a a c o 
a © c o 
ooc© 

C n t- © 
N* 43 I— Q 

4) m c © 
O O o Q 
in © r- 45 


o c 
o c 
oo 
a a 
o o 
cm a 
45 N“ 

rH IT 


©ncoooo© 


«-* TO CO 

KoE-po go -P 
P-* Ph Eh E-t Eh Q. Cl 



TABLE III.. - A SUMMARY OF PVT RELATIVE ERRORS; FROM FIGURES 1 
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A comparison of the data of Vennix to Vennix's equation of state. CS-60063 

e* 2 

The Strobridge equation has been overextended to the range up to 1000 K. and 100 MN/m . 

c** 

Equation should not be used below 77 K. 
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TABLE IV. - OPERATIONS SHEET FOR SUBROUTINE GASP 


COMMON/PROPTY/KU, DL, DV, HL, HV, S, SL, SV, CV, CVL, CW, CP, CPL, CPV, GAMMA, GAMMAL, GAMMAV, C, 
CL, CVP, MU, MUL, MUV, K, KL, KV, SIGMA, EXCL, EXCV, EXCESK 
REAL MU, MUL, MUV, K, KL, KV 

CALL SETUP(NAMGAS) Name the fluid: HE for Helium, CH4 for methane, NE for Neon, N2 for Nitrogen, CO for 

Carbon Monoxide, 02 for Oxygen, AR for Argon, C02 for Carbon Dioxide. 

CALL GASP(K3, KP, T, P, D, H, KR) 


\ 


1 I 

i i 


i i 


\ 


\W 

V \ x 


\ 


\ 


x \ \\ 

ww 

\ 

\ 


1 


L 


Region 

KR=0 Unknown, check KR returned 
KR=1 Saturation 
KR=2 Liquid 

, KR=3 Gas and or fluid 
Enthalpy, J/g 
\ Density, g/cm 3 
\y- Pressure, MN/m 2 
^Temperature, K 


\ u Thermodynamic and Transport properties* 



\ 


KP=0 Only P, p, T returned 
KPfcl H Enthalpy, J/g, (H), (HL), (HV) 

RP*2 S Entropy, J/(g)(K), (S), (SL), (SV) 

KPt=4 Cv Specific heat at constant volume, J/(g)(K), (CV), (CVL), (CW) 

Cp Specific heat at constant pressure, J/(g)(K), (CP), (CPL), (CPV) 
y Ratio of specific heats, Cp/Cv, (GAMMA), (GAMMAL), (GAMMAV) 
c Sonic velocity, cm/sec, (C), (CL), (CVP) 

KP=8 p Dynamic viscosity, g/(cm)(sec), (MU), (MUL), (MUV) 

KP=16 k Thermal conductivity, W/(cm)(K), (K), (KL), (KV) 

KP=32 a Surface tension, dyne/cm, (SIGMA) 


- Input specification of independent properties 


KS=1 

KS=2 

KS=3 

KS=4 

KS=5 


p = f(T,P);*fWS$fl 1T,P fin$ p 
P=f(T,p); giWn T,p find P 
T=f(p ,p)j given P,p find T 
T,p = f(P,H)j_ given P,H find T,p 
T, p « f(P,S); given P, S find T,p 
Reset KR/1 each time GASP is called to be assured of nonsaturation calculations. At the beginning of the user's program, 
set KU=1 for units indicated above. For other units, see Table V. 

SAMPLE PROBLEM 

COMMON/PROPTY/ i(as indicated above) 

DATA NAMGAS/2HN2/ 

CALL SETUP(NAMGAS) 

KU=1 

CALL GASP(2, 31,80, P, 1.254, H,KR) 

The program has been asked to find for nitrogen the pressure P, corresponding to a density of 1. 254 g/cm 3 and a tempera- 
ture of 80 K. Furthermore, the program will return values for H, S, Cv, Cp, y, C, p, and k. The values of H and P are 
returned through the call vector and the remaining values are returned through the COMMON statement. 

* 

For example: KP= 1+2=3 will return enthalpy and entropy; KP= 1+8+16=25 will return enthalpy, viscosity, and thermal 
conductivity; KP= 1+2+4+8+16=31 returns everything except cr; KP=2 will return entropy. 



TABLE V. - UNITS SPECIFICATION KU 


Physical quantity 

KU=1 

KU=2 

KU=3 

Temperature 

K 

K 

R 

Density 

gm/ cc 

gm/cc 

lb m/ft J 

Pressure 

MN/m 2 

atmospheres 

psia 

Enthalpy 

joule /gm 

j oule/ gm 

Btu/lb m 

Entropy, specific heat 

joule /gm-K 

joule/gm-K 

Btu/lb m-R 

Sonic velocity 

cm/ sec 

cm/sec 

ft/sec 

Dynamic viscosity 

gm/ cm-sec 

gm/cm-sec 

lb m/ft-sec 

Thermal conductivity 

j oule/ cm-sec-K 

joule/cm-sec-K ' 

Btu/ft-sec-R j 

Surface tension 

: dyne /cm 

dyne /cm 

ibf/ft j 

L — i 


KU=4,5 permit the user to work in other units; however, the proper conver- 
sions must be entered into BLOCK DATA. 
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TABLE VI. - PROGRAM ASSEMBLY 

1» Set up the return of thermodynamic and transport properties 

COMMON PRPTY/KU ,DL,DV,HL,HV,S,SL,SV,CV,CVL,CW,CP, CPL , CPV , GAMMA , 

GAMMAL , GAMMAV , C , CL , CVP ,MU ,MUL ,MUV,K,KL,KV, SIGMA, EXCL ,EXCV ,EXCESK , 
REAL MU,MUL,MUV,K,KL,KV 

2. Specify the name of the fluid for which properties are required; e.g., 

DATA NAMGAS / 3HCH4/ . 

3. Call SETUP (NAMGAS). This overlays the proper coefficients. See 

table II. 

4. Establish the units you are working in KU= (specify 1, 2, or 3). See 

table VI. 

5. Add your program here. 

6. Add GASP package or SUBROUTINE SETUP (NAMGAS), BLOCK DATA, and other 

required subroutines. 

Data 

GASP Package or SUBROUTINE 
SETUP, BLOCK DATA, and 
other subroutines 

Main Program + Subroutines 

KU=1, 2, or 3 

DATA NAMGAS / 3HCH4 / 

CALL SETUP (NAMGAS) 

COMMON PRPTY, etc. 

Program Dimensions, Common, 
etc. 

Control cards 
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Figure 1. - Relative error in density for the methane data of Vennix, refer- 
ence 57, as a function of reduced density; calculated values from Vennix's 
equation of state, references?. cs _ 60065 



Figure 2, - Relative error in pressure for the methane data of Vennix, refer- 
ence 57, as a function of reduced pressure* calculated values from Vennix's 
equation of state, reference 57. cs-6oo66 
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Figure 7. - Relative error in density for the methane data of Jansooue et. al f 
reference 56 as a function of reduced density; calculated values from 
Bender's equation of state, reference 22. cs-60071 
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Figure 8. - Relative error in pressure for the methane data of jansooue et. al, 
reference 56, as a function of reduced pressure; calculated values from 
Bender's equation of state, reference 22 . cs-60072 
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Figure 10. - Relative error in density for the oxygen data of Weber, refer- 
ence 31, as a function of reduced density; calculated values from Benders 
equation of state, reference 22. cs- soots 



Figure 11. - Relative error in pressure for the oxygen data of Weber, 
reference 31, as a function of reduced pressure,- calculated values 
from Bender’s equation of state, reference 22. cs- 6 oo 76 
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Figure 12. - Relative error in temperature for the oxygen data of Weber, 
reference 31, as a function of reduced temperature; calculated values 
from Bender's equation of state, reference 22. ..cs-60077 



Figure 13. - Behavior of the specific heat at constant volume, 
Cv, calculated using d 2 P/dT 2 relative to C v 
computed from NTWO using (AU/AT)p. Pressure * 15 atm. 
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Figure 14. - Thermodynamic and transport properties for fluid methane, cs -60059 
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Figure 15.- Selected output from GASP for various isobars of fluids oxygen 
and helium (0. 101 325 * ATM. * MN/M^). 
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